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ABSTRACT 

The theoretical predictions of big bang nucleosynthesis (BBN) are dominated by 
uncertainties in the input nuclear reaction cross sections. In this paper, we examine the 
impact on BBN of the recent compilation of nuclear data and thermonuclear reactions 
rates by the NACRE collaboration. We confirm that the adopted rates do not make 
large overall changes in central values of predictions, but do affect the magnitude of the 
uncertainties in these predictions. Therefore, we then examine in detail the uncertain- 
ties in the individual reaction rates considered by NACRE. When the error estimates 
by NACRE are treated as la limits, the resulting BBN error budget is similar to those 
of previous tabulations. We propose two new procedures for deriving reaction rate un- 
certainties from the nuclear data: one which sets lower limits to the error, and one 
which we believe is a reasonable description of the present error budget. We propagate 
these uncertainty estimates through the BBN code, and find that when the nuclear 
data errors are described most accurately, the resulting light element uncertainties are 
notably smaller than in some previous tabulations, but larger than others. Using these 
results, we derive limits on the cosmic baryon-to-photon ratio r/, and compare this to in- 
dependent limits on r] from recent balloon-borne measurements of the cosmic microwave 
background radiation (CMB). We discuss means to improve the BBN results via key 
nuclear reaction measurements and light element observations. 
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1. Introduction 

Big bang nucleosynthesis (BBN) has long played a key role in the standard cosmology. The 
concordance between the predicted and observed light element abundances (Walker et al. 1991; 
Olive, Steigman, &; Walker 2000) gives confidence that the basic framework is working down to 
epochs of t ~ 1 sec, and z ~ lO^'^. The theory-observation comparison also probes particle physics 
in the early universe (Sarkar 1996; Olive & Thomas 1999; Lisi, Sarkar, &: Villantc 1999), and 
quantifies the allowed range in the baryon-to-photon ratio r] = n-Q/n^ (Fields &: Olive 1996, Fields 
et al. 1996, Schramm & Turner 1998) and thus the cosmic baryon density parameter ^l-^h? (where 
= pj/pcrit = 8TrGpi/3H^, and Hq = lOO/i km s~^ Mpc~^). 

Over the past decade, a major thrust of research in BBN has been towards increasing the rigor 
of the analysis. On the theory side, the key innovation was to calculate the errors in the light 
element predictions in a systematic and statistically careful way. This was done using Monte Carlo 
analyses (Krauss & Romanelli 1990; Smith, Kawano, & Malaney 1993; Kernan &; Krauss 1995; 
Hata et al. 1995; Fiorentini et al. 1998; Nollett & Buries 2000), which account for nuclear reaction 
uncertainties and their propagation into uncertainties in the light element abundance predictions. 
These calculations are essential because they allow for a careful statistical comparison of BBN 
theory with observational constraints; in addition, they point the way toward improvements in the 
theory calculation. Just as BBN theoretical and observational ingredients have been sharpened, 
so have the tools for comparing the two. Statistical tests involving likelihood analyses (Hata et 
al. 1995; Fields, Kainulainen, Olive, & Thomas 1996; Lisi, Sarkar, Sz Villante 1999) have been 
developed and performed. 

A focus on rigor in BBN will soon be rewarded, as cosmology moves toward a precision era. 
Specifically, the observations of the anisotropics in the cosmic microwave background (CMB) will 
allow for very precise determination of cosmological parameters, including Q^h^ (or, equivalently, 
7]). A comparison of rj as determined by BBN and the CMB will provide a fundamental test of 
cosmology. In addition, as pointed out by Schramm &; Turner (1998) and as shown in detail below, 
increasingly accurate CMB data can ultimately transform BBN into a much sharper probe of the 
early universe and of chemical evolution. 

Along with the neutron lifetime, eleven key nuclear reaction rates represent the dominant 
sources of error the BBN calculation (Smith, Kawano, & Malaney (1993), hereafter SKM; and §4 
below). Thus, the choice of thermonuclear reaction rates, and their uncertainties, will determine 
the accuracy of the final predictions. Most recent work on BBN has used the Caughlan &; Fowler 
(1988) compilation, with updates due to SKM and others. The SKM error budget has been the 
standard for all of the subsequent Monte Carlo work except for that of Nollett &; Buries (2000), 
who create their own rates from the nuclear data, but do not present the thermonuclear rates by 
themselves. The recent work by NACRE (Angulo et al. (1999))^ represents a significant new effort 



^http : / /pntpm . ulb . ac . be/nacre . htm 
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to critically evaluate nuclear cross section data and derive thermonuclear rates. Moreover, the 
NACRE compilation includes estimates of uncertainties in the reaction rates. 

We therefore have examined the impact of the NACRE rates on BBN. This issue has also 
been studied by Vangioni-Flam, Coc, & Cassc (2000). These authors considered the impact of the 
error range estimated by NACRE by calculating the effect of individual rates on the primordial 
abundances. Vangioni-Flam, Coc, &: Casse (2000) then estimated the effect of the ensemble of 
rates by placing all rates at their maximum and minimum variation. Here, we will extend this 
work in several ways. We will use a full Monte-Carlo treatment, which allows for a quantitative 
statement about the propagation of the errors. In addition, we examine the impact of different 
error assignments: (1) those of the NACRE group themselves, which give useful estimates of the 
range of uncertainty but are not defined in a uniform manner; as well as (2) the results of our own 
error analyses, which we derive from the nuclear data according to simple but uniform procedures. 

This paper is organized as follows. We will describe the NACRE recommended rate uncertain- 
ties in §2. In §3 we compare the NACRE uncertainty estimations with other estimates based on the 
nuclear reaction data. These different error estimates are used to predict light element uncertainties 
via Monte Carlo BBN calculations (§4). The results are compared with observations via likelihood 
analyses in §5, and discussed in terms of the predictions for rj. We discuss the agreement of our 
predictions with current CMB data in §6, and anticipate the impact of high-quality data which will 
arise from future CMB space-based experiments. Our conclusions are summarized in §7. 



In an effort to update the thermonuclear reaction rate compilation by Caughlan &: Fowler 
(1988), the NACRE collaboration has presented a detailed analysis of 86 charge induced nuclear 
reactions. Out of these reactions we will discuss the 7 reactions NACRE has in common with the 11 
reactions from SKM(1993), whose errors dominate the uncertainties in the abundances. (see §2.3) 



The nuclear reaction inputs to BBN take the form of thermal rates. These rates are computed 
by averaging nuclear reaction cross sections over a Maxwell-Boltzmann distribution of energies. 
The thermonuclear reaction rate at some temperature T, is given by the following integral: 



where A^^ is Avogadro's number, v the relative velocity, gives the reduced mass of the nuclei, ks 
Boltzmann's constant and cr{E) yields the cross section at center of mass energy E. 



2. The NACRE Compilation 



2.1. 



Reaction Rate Formalism 




(1) 
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The charge induced cross sections can be decomposed into 

E 

where S{E) gives the astrophysical S-factor, C is the Sommerfeld parameter, defined by 

C = z,zJ|!)n=z,zJ|!ViSV. (3) 



hcj \vJ \hc J \ 2E J 

Here the ZiS are the charges of the nuclei in units of the proton charge e. 

For neutron induced reactions, the cross section can be written as fohows: 

a{E) = (4) 
where R{E), the i?-factor, can be a slowly varying function of energy, and is similar to an S-factor. 



2.2. Recommended Thermonuclear Rates 



Generally speaking, NACRE follows the standard path (outlined in the previous section) for 
going from nuclear data to thermonuclear rates. NACRE has compiled an extensive tabulation of 
nuclear reaction cross section data, which is available in a very convenient form online at the NACRE 
website http://pntpm.ulb.ac.be/nacre.htm. For each experiment, the NACRE collaboration 
has organized cross section data into useful tables, providing information on the measured cross 
sections and their total error as a function of energy. These tabulations axe invaluable for the 
analysis of error propagation. 

It is important to have a consistent method of analyzing cross sections and computing reaction 
rates. The data sets provided by NACRE over the relevant energy ranges for BBN have large 
variations in both quality and quantity. Thus, NACRE critically evaluates the experimental data 
and gives different data sets different weights in determining rates. The NACRE collaboration 
considers only charged induced reactions, thus leaving out several important reactions for BBN. 
Therefore we have to include the neutron-induced reactions and one deuteron-induced reaction by 
the SKM compilation. The reactions used from each are shown in Table 1. For completeness we 
will describe the approach for the neutron-induced reactions as well. 

Once the cross section data (and errors) are gathered, they are put into S'-factor or i?-factor 
form. NACRE fits these factors to analytic functions. For non-resonant data, a polynomial of order 
2 to 3, is used. Some data are fit from theory (e.g., NACRE uses Kajino's (1986) fits for t{a, ^Yhi 
and '^Hc(a, 7)^Be), while for resonant data a Brief- Wigner or R-matrix fit is performed (e.g., for 
t{d, n)'*Hc). The fits are required to be a good representation of the data by a analysis, taking 
precedence over agreement with theory. 
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The evaluated cross section fits, S'fit(-E') are used to determine thermal rates (via (1)) for a 

grid of temperatures, T, ranging from Tg = .001 — 10, where T = Tg x lO^K. These results are 
subsequently fit to a prescribed function of Tg. One should note that the NACRE fits do not rcficct 
the analytic approximation for the thermal rates used by Caughlan & Fowler (1988). Namely, 
Caughlan &; Fowler put non-resonant thermal rates in the form (^^- CjT^^'^) cxp(a/Tg^'''^). NACRE 
modifies this approach, keeping the same exponential dependence, but changes the prefactor from a 
polynomial in Tg^^ to one in Tg: c'jT^) exp{a' /Tg^^). The main reason for the form of their fit 
is to get fast convergence to the numerical data. In some cases (e.g. '^He(d, n)^He and '^Li(p, Q;)^He) 
additional factors are used to improve the fit to the numerical results. 

Table 1: Key Nuclear Reacti ons for BBN 
Source Reactions 
NACRE d{p,-rflle 

d{d,p)t 
t{d, n)^He 

3He(Q!,7)^Be 
^Li(p, «)4He 

SKM p{nrf)d 

^Re{d,pf}ie 

'^Be{n,pyLi 
This work ^He(n,p)i 

PDG Tn 



As noted above, some of the rates are not provided by NACRE. In these cases, the SKM rates 
as indicated in Tabic 1 are used. One of these, '''Be(n,p)''Li, is a n-capture reaction for which a 
large amount of data is available. The deuteron- induced reaction (^He(d, p)^Hc), is fit as a charged 
particle reaction using the Caughlan &: Fowler prescription, as discussed in the previous paragraph. 

Several reactions deserve special mention. As noted by SKM and emphasized recently by 
Nollett & Buries (2000), the p{n,^)d reaction suffers from a lack of data in the BBN energy 
range. Also, p{n,j)d has only 4 data points (not available when SKM did their study) in the 
relevant energy range < 1 MeV. Fortunately, this reaction is well-described theoretically. Here we 
follow both SKM and Nollett &: Buries, by adopting the theoretical cross sections of Hale ct al. 
(1991), which provide an excellent fit to the four available data points by Suzuki (1995) and Nagai 
(1997). Nevertheless, despite the present agreement between theory and data, the importance of 
this reaction-which controls the onset of nucleosynthesis-demands that the theoretical cross section 
fit be further tested by accurate experiment. We urge further investigation of this reaction. 

Since SKM, Brune et al. (1999) have added new and very precise data for ^He(n,p)t (see Figure 
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la).^ This has greatly reduced the uncertainty in this reaction. In order to use these data, we have 
refit the R factor in the manner of SKM and Brune et al., using a third order polynomial in v and 
the entire world data set. We arrive at fit parameters very similar to those of SKM and Brune et 
al., with Na{(7v) = 7.3546 x 10^(1 - 0.7757TJ/^ + 0.5376r9 - O.lOlSTg^^) cm^ g-^ 

The weak reactions which govern n <^ p interconversion also deserve mention. These reactions 
are a strong function of temperature, but the rates can be scaled to a single laboratory measurement, 
the neutron lifetime r„. Thus, the uncertainty in the weak rates is set by the uncertainty in r„. 
Here, we have adopted the value recommended by the Particle Data Group, r„ = 885.7 it 0.8 s. 
This value reflects the recent and very precise measurement of Arzumanov et al. (2000), and has 
considerably reduced the errors of the Particle Data Group world average from the old determination 
of 886.7 ± 1.9 s (Groom et al. (2000)). 

It is also important to include the small but non-negligible corrections to the tree-level weak 
rates. These include a number of contributions, such as radiative corrections (Dicus et al. (1982); 
Heckler (1994); Esposito, Mangano, Miele, &: Pisanti (2000a)), nonequilibrium neutrino heating 
during 6=*= annihilation (Dodelson Sz Turner (1992)), and finite nucleon mass and finite temperature 
effects (e.g., Seckel (1993); Kernan (1993); Lopez Sz Turner (1999); Esposito, Mangano, Miele, & 
Pisanti (2000a)). These corrections have been implemented in our code as has been described in 
detail in Olive, Steigman, & Walker (2000). As discussed there, the differences in the ^He yields 
between our code and that of Lopez & Turner (1999) are 0.0001 ± 0.0001 for 1 < r?io < 10. 

2.3. Uncertainty Limits 

NACRE uses data uncertainties in two ways, for the evaluation of mean rates and the evaluation 
of errors. For the mean ("adopted") rates, they use the data errors to weight sets via a analysis. 
They also use the data errors to estimate the uncertainty range for these "adopted" rates. This 
is not done in a strict statistical way, and thus is not presented as, e.g., "Icr" error, but rather as 
"high" and "low" limits to the thermal rates. 

The "high" and "low" limits are derived in different ways for different reactions. In the case 
where there happens to be two discrepant data sets, NACRE performs a analysis for each of 

the two data sets, adopting the set with the larger S'-factor as their "high" limit, and the set with 
smaller S'-factor as their "low" limit. The "adopted" value for their mean is simply the unweighted 
average of the "high" and "low" limits (e.g., d(p, 7)^He). 



^Note that in all figures having logarithmic vertical scales, errors have been properly propagated to reflect the log 
nature of the plot. 
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3. A Comparison of Nuclear Reaction Uncertainty Estimations 

The power of BBN theory to test and constrain cosmology and particle physics derives from 
the ability to calculate accurately the mean values and uncertainties in light element abundances 
as a function of rj. The Monte Carlo analyses of BBN which provide these results are in turn only 
as good as the input nuclear and weak error budget. It is thus crucial to make the most accurate 
estimates possible, in a way that is appropriate for the cross section data sets and the thermal rate 
compilation one uses. 

Given the importance of the error propagation, we will explore different methods of doing this 
for the NACRE compilation. We first examine NACRE's own, "high/low" limits. As these are 
not uniformly derived, we also consider two other error budgets, which illustrate issues in error 
propagation and subsequent impact on BBN. In this section, we will describe the error propagation 
methods, and our implementation of them. We will then examine the BBN results using these 
errors in the following section (§4). 

3.1. Tailoring Rates and Uncertainties for BBN 

Since the error analysis of the rates by NACRE is not done in a strict statistical way, it 
becomes important to develop a general, statistically sound method for determining an accurate 
mean and error for the application to BBN. To do this, we will compare the NACRE theoretical 
S*- factor fits to the data sets they use. For the theoretical fits Sfn (E) , we for the most part use the 
analytical forms published in the NACRE paper. When these were unpublished, NACRE kindly 
made available their fits (or in the case of n)^He, a tabular form).^ These fits are valid for 
different energy ranges depending on the reaction, but in all cases cover an energy band which 
includes and extends beyond the energies needed for BBN. 

Our goal is to assess the uncertainties in the cross sections needed for BBN. To do this, we will 
use a analysis to describe the goodness of fit of the theory to the data, given the experimental 
errors. However, it turns out that the NACRE cross section fits do not give the minimum over 
the energy range for which the fits are valid. As we will see the shape of the fits represent the data 
well, but the normalizations are not always precisely those which minimize x^- This reflects the 
fact that the NACRE S-factors are in some cases designed to fit the cross section data over larger 
energy ranges, and in some cases were chosen to be a compromise between conflicting sets of data. 
To hone the NACRE rates for use in BBN, we allow for a shift in the overall normalization of the 
S'-factor fits. This allows for a simple means of refining the rates for BBN, while maintaining the 
shape of the cross section fits and thus the shape of the NACRE thermal rates. Since the fits are 
valid for energies beyond the BBN range, their shapes are more strongly constrained than they 



^We are particularly indebted to C. Angulo and P. Descouvrement (2000, private communication) for their help. 
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would be if we considered only the BBN energies. 

Our procedure is as follows. For a given reaction, we have a set of measured S-factors 
{Sohsi^i)} nieasured at energies {E^}, where i e 1, N. The S-factors have errors {cTj}, while the 
energies are measured with negligible error. The measured S'-factors are described by theoretical 
fits S{it{E), which combine theory and phenomenology as discussed above. 



Table 2: Reduced for Key BBN Reactions 



Reaction 


a-1 


# Points 


Amin 




d{p, 7)^He 


-0.076 


72 


987 


13.9 


t{d,n)^lle 


0.000 


203 


295 


1.46 


d{d, nfRe 


+0.018 


52 


312 


6.11 


d{d,p)t 


+0.006 


92 


425 


4.67 


t{a, j^'Li 


+0.058 


30 


8.42 


0.290 


=^He(a,7)^Be 


-0.067 


131 


242 


1.86 


^Li(p, a)^He 


-0.045 


126 


374 


2.99 


^He(n,p)t 


0.000 


165 


56.9 


0.347 


3He(d,p)^He 


+0.0694 


73 


171 


2.34 


p{n,^)d 


-0.012 


4 


0.901 


0.300 


'^Be(n,p)^Li 


+0.010 


65 


38.9 


0.608 



We allow for a different normalization a between the "true" theory curve and the fit curve. 
That is, we write the theory curve as 5th (-E) = ctSfit{E), where we expect a to be near, but not 
necessarily equal to, unity. Taking the data errors to be Gaussian, the data and theory agreement 
is quantified by the value of chi-squared: 

X =x{a) = 2^ ^2 = ^2 (5) 



The renormalization is determined by minimizing %^ with respect to a, giving the best-fit value of 
a to be 

which will be close to unity if the experimental data scatters evenly around the fit curve, as expected. 



The renormalizations appear in Table 2. We see that in general, the shifts are small, less than 
10% for all cases. The effect of renormalization can be seen graphically by comparing the solid 
and long-dashed curves in Figure 1. As the figure illustrates, the changes are slight-by eye, the 
goodness of fit of the original NACRE fits is similar to that of the renormalized fits. However, the 
improvement of the fits is significant enough to justify the addition of a parameter (a). We will 
see in the next section (§4.1) that the renormalizations have a mild effect on central values of the 
BBN predictions, but the shifts are nonetheless significant since they are comparable to the level 
of the uncertainties. Thus, we will show the effect on BBN of adopting the renormalizations. 
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3.2. NACRE High/Low 

As discussed above (§2.3), for each reaction, NACRE presents an "adopted" thermonuclear 
rate A, as well as a "high" rate Ah and "low" rate A^. Each of these is of course a (strong) function 
of T. NACRE describe the high/low rates as "lower and upper limits" to the adopted rates, but 
do not present them as statistically defined limits, such as la or 2a ranges. Thus, if we wish 
to use these limits, we must first determine what statistical weight to give them, and then set a 
prescription for using these limits in the Monte Carlo procedure. 

We thus turn to the question of assigning statistical meaning of the NACRE high/low limits. 
The uncertainties are conveniently quantified by the fractional or relative error 5i = |Aj/A — 1|, 
with i = high and low In principle, for a given reaction rate the high and low errors need not be 
equal, but in practice we find that the high and low fractional errors are nearly identical (certainly 
within our desired accuracy) over the entire temperature range. Thus, we can write the high/low 
limits as Aj(T) = A(T) [lib(5(T)], where 5 = 5±. This simplifies the analysis slightly. Figure 2 plots 
the fractional errors 6{T) for the seven NACRE reactions important for BBN. For comparison, we 
have also plotted fractional errors given in the SKM compilation.^ Several trends are apparent 
in the NACRE curves. We see that in general, the NACRE high/low fractional errors vary with 
temperature; this is in contrast with the SKM rates, for which only 2 of the 12 rates are assigned 
a temperature-dependent fractional error. 

Comparing the high/low limits with the la SKM errors, we see that the errors are of the 
same order of magnitude. There is no clear trend as to which compilation has the larger errors. 
For example, the NACRE errors for i(a,7)^Li are lower than the SKM rates (due in part to the 
addition of new, accurate data by Brunc ct al. (1994)), while for the mirror reaction ^He(a, 7)'''Be, 
the NACRE limits are larger than those of SKM. For the purpose of comparison with our more 
rigorous results which follow, we have taken the simple approach of assigning the NACRE high/low 
limits to be la errors. Given the variations in the NACRE rates and their comparison with the 
SKM variations, this approach is certainly qualitatively correct. A similar approach was taken by 
Vangioni-Flam, Coc, & Casse (2000), who examined the impact of the NACRE rates and ran cases 
using the high/low limits to compare with the central rates. 

While the NACRE assignment of limits to their rates invites quantification in a full Monte 

Carlo code, the somewhat arbitrary nature of the specific assignment of la errors to the limits leaves 
question as to the robustness of the results. While here we have assumed a Gaussian distribution 
for these assignments, other choices such as a log-normal distribution are also possible (Vangioni- 
Flam, Coc, & Casse (2000)). Thus, we now turn to other error assignments which are derived from 
the underlying cross section data. 



"^In the BBN code, for the three "SKM-only" strong reactions (Table 1) as well as for ^He(n,p)t, we apply the same 
analyses as in the case of the NACRE rates, with the exception that NACRE's high/low error estimate is replaced 
by the SKM recommendation. 
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3.3. Minimal Uncertainties: Ax^ = 1 

To derive errors for NACRE rates in a more statistically well-defined way, we now perform our 
own analysis of data and propagation of errors. In this section, we will treat all cross sections errors, 
from all data sets, as independent. We thus ignore any correlations between the measurements, 
even within a given experimental run. Clearly, this is not realistic, but this procedure has the 
virtue of simplicity and will give the limiting case of the smallest possible error achievable with 
nuclear data. This will also serve as a useful starting point for a more realistic analysis in the next 
subsection. 

We will compute the variance in the S'-factor in the form of a fractional error. This is also an 
approximation, but simplifies the analysis and the limit that results. Physically, this corresponds 
to the assumption that the large amount of data determines the shape of the curves quite well, so 
that the dominant uncertainty is not in the shape but in the overall normalization. The data sets 
used are the same as NACRE; the online versions include data and errors which is greatly helpful 
in this process. 

We can determine the fractional error in the theory fit by varying the normalization until 
= Xmin + which occurs when 

ai = ai: — , (7) 

Thus, we have a la variation when a is changed from its minimizing (renormalized) value by a 
fractional error 

/ = ai/a - 1 (8) 

The theory-data fit thus has a it la error ath = /<5th = afSm- 

We have found these "Ax^ = 1" errors for the 11 strong reactions; these appear in Table 3. 
In the Table, we see that the errors are very small, < 1%. This can be seen to follows from the 
large number of data points. Of course, we have in this case we have assumed each measurement 
in each experiment to be independent of all others. We have thus ignored the correlations among 
the errors. These problems are highlighted when we examine the questions of goodness of fit. 

The goodness of fit is quantified by xi — X^/^j the x^ pcr degree of freedom with = — 1. 
Table 2 gives the reduced x^ for the 11 key strong reactions. For normally distributed and correctly 
estimated errors, one expects Xu to lie close to 1. In fact, this is not what is found for any of 
the reactions. In four cases, the reduced x^ is significantly less than 1, while in seven cases x^ 
significantly exceeds 1. In Figure Id, we see the data and fit for t{a, 7)^Li, an example of a reaction 
with x^ < 1- In this case, we see that all data are within la of the fit, and indeed some lie 
essentially on the fit, leading to the small reduced x^ and suggesting that the experimental errors 
may be overly conservative in this case. 

For the cases in which x^ > Ij we follow the practice of the Particle Data Group and we 
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increase the fractional error by an amount \/x^. This scaUng encodes the presumption that at 

least one of the experiments has underestimated its uncertainties. Where the data are discrepant, 
this is most assuredly the case. For the cases in which the xi < li no correction is made. 

Let us examine more closely the large Xu cases, for which the goodness of fit is poor. A 
particularly egregious case is d{p, 7)"^Hc, with xt — 13.9, which appears in Figure lb. The reason 
for the poor fit is apparent: the data for this reaction are internally inconsistent in the 0.01 to 0.05 
MeV range. Specifically, the recent measurements of Schmid et al. (1995) do follow the fit trend, 
but the older and less precise data of Griffiths, Lai, &; Scarfe (1963) lie above these data, outside 
of the range of the quoted errors. In our fitting procedure, the small errors in the Schmid et al. 
(1995) data dominates the fit normalization, but the systematic offset of the Griffiths, Lai, &; Scarfe 
(1963) data leads to a poor x^- 

Another case with a high Xu '^He(Q!, 7)^Be, which appears in Figure Ic. This reaction has 
been wcll-studicd in the context of solar models, where it controls the intensity of '^Bc neutrinos 
and thus is of great importance for the solar neutrino problem. Studies of this reaction are reviewed 
by Adelberger et al. (1998), who note the discrepancy between 7-capture and ^Be activity mea- 
surements. Although the discrepancies between the data sets are not as glaring as in the d{p, 7)^He 
case, the systematic offset between data sets is nevertheless clear. 

The source of the systematic error traces back to the difficulty in making an absolute cross 
section measurement. While a given experiment can accurately determine the relative values be- 
tween different energies, measuring the absolute cross section values involves calibration of beam 
intensities, beam energy losses in the target, detection efficiencies, solid angle uncertainties, etc. 
Thus, the absolute cross sections from each experiment carry a normalization error, which is com- 
mon to all points from the experiment. As emphasized by Buries, Nollett, Truran, & Turner (1999) 
and Nollett & Buries (2000), these normalization errors in the cross sections affect all of the BBN 
reactions, and are properly treated as correlated errors. 

Thus, our assumption of independent errors, and the Ax^ = 1 method, while simple and illus- 
trative, has greatly underestimated the true error budget, in which correlations play an important 
and even dominant role. We now turn to a method which includes these effects and treats the data 
in a more uniform manner. 

3.4. Sample Variance 

We now seek a new error estimator which accounts for correlated errors due to absolute nor- 
malization uncertainties in the cross section measurements. These dominate the uncertainties for 
many reactions. Furthermore, these uncertainties are not reduced by taking more data within a 
given experiment; rather, each experiment effectively represents one attempt to measure the ab- 
solute normalization, independent of the number of points measured in the experiment. Thus, 
the normalization errors for a given reaction are only lowered by combining many independent 
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experiments. Consequently, the overall error budget in the cross sections, and thus in the inferred 
thermonuclear rates, is much larger. 

We note that the systematic errors in absolute cross section normalization are quoted as 
percentages. Thus, the use of fractional errors is justified when these errors dominate. We will thus 
cast our uncertainty determinations in the form of a fractional or normalization uncertainty (as in 
the previous section), but now with more justification. 

In constructing an error estimate which accounts for correlated normalization errors, the data 
sets with discrepant normalizations will provide a useful test case. Our goal is to find a well-defined 
procedure which will give errors that automatically account for systematic offsets among data 
sets, leading to an error range which accommodates both sets. We therefore adopt a (fractional) 
error determined by the weighted sample variance, as follows. The notation is the same as that 
introduced in the previous section. In particular, we again adopt a theory normalization a given 
by 5th = ct'S'fitj and use the data to determine the value of a and its error. The mean value a of 
the normalization is determined again by minimization. Given this, we can find the fractional 
difference 

r _ 'S'obs(-E^i) - StYijEj) _ So\,s{Ei) _ , . 

' ~ SUE^i) ~ aS^t{Ei) ^ ' 

between each data point and the theory, and the fractional error 

ei = (Ti/Sti,{Ei) = ai/aS^t {Ei) (10) 

in each data point. Given these quantities, the theory normalization a which minimizes i^l- 5) 
is also the value of a which guarantees that the weighted average theory-data difference satisfies 
((5) = Yl^i/^l = 0- The error, as measured by the weighted sample variance aa, is now 
determined as the weighted average of Sf: 



2 _ ^1/^ _ J:^ [SoUE^) - S,,,{E,f lal 



Note that if a typical data point S'obs hcs o away from the theory point iSth, then (T„ ~ (j/5th, the 
fractional error. However, since the cross sections are a function of energy and are not constant, 
the weighting will favor points which have a small fractional error ailSxh{Ei)- 

We now turn to the behavior of eq. (11) in the case of discrepant data sets. For a simple but 
illustrative example, consider a reaction for which there are two data sets which are inconsistent 
with each other. Let each data set contain TV points, with comparable errors. Let the first set have 
measurements with So\,^{Ei) = Q!i5fit(-E'j) -|- 5j for i G (1, . . . , A?^), so that the data scatters about 
the fit curve with a normalization discrepancy (systematic error) ai, and the errors which scatter 
as {bi) = and {bf) = af. Similarly, the second data set has points with Sohs{E,i) = a2S{it{Ei) +Si 
for i G (A^-|- 1, . . . , where a2 / ai indicates a systematic discrepancy, and {6f) = af. Finally, 

as described above, we allow for a theory normalization error: iSth = a-S'fif Then the expected (x^) 
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is minimized for the normalization 

i.e., a weighted average of the two factors Uj. This result is of course quite reasonable. In this case, 
we have 

{x'{a)) = 2N + {Aa/2fY.^^^ (13) 

1=1 '^i 

We see that differs from the number v = 2N — 1 ~ 2N of degrees of freedom by an amount 
which depends on Aa, the difference in the systematic error. In this case, therefore, the degree to 
which exceeds u, or xt exceeds 1, is just a measure of the systematic discrepancy between the 
data sets. Again, this is as expected. 

For the special case when S'gt = 5*0 and crj = (Tq, both constants independent of energy, eq. 
(12) reduces to just a = {ai + a2)/2, the simple average. Furthermore, in this case the weighted 
sample variance of eq. (11) takes the simple form 

al = {ao/Sof + {Aa/2f (14) 

where Aa = q;2 — ai. Again, this result is quite reasonable: the fractional error in the normalization 
is the quadrature sum of a typical fractional measurement error (Tq/Sq and the systematic error 
as estimated by (a2 — q;i)/2, the difference in normalizations. Thus, when the systematic error 
dominates, then ao/So <C Aa we have an overall fractional error cTq — Aa/2, so that the range 
a = a lb (Tq, corresponds precisely to the range ai < a < a2. This is precisely what one expects 
if the systematic error dominates; our expression automatically gives an answer which spans the 
range between the discordant data sets. On the other hand, if the systematic error is negligible, 
then we have uq/Sq S> Aa, and (Tq, ~ (Jo/So, a typical measurement error. In this limit, then, the 
fractional error is such that it encloses the la error band about the (normalized) fit curve. 

Table 3 shows ctq for the 11 strong BBN reactions. We see that these results are either 
comparable to or smaller than the SKM results. The largest uncertainties appear when incompatible 
data sets are present. As expected, the systematic offsets are accounted for by an increase in the 
errors. 



3.5. Other BBN Error Studies 

Other studies have been made which determine errors in the nuclear inputs for BBN and 
propagate them to derive uncertainties in BBN abundance predictions. The SKM compilation 
of nuclear uncertainties has been very widely used for BBN Monte Carlo studies. This work 
was extremely influential, as it was the first to catalog the most important BBN reactions, to 
systematically collect, and display the available cross section data for them, and to derive a set of 
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Table 3: RuacLioii Errors 
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0.07 
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0.00621 
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0.09 



errors for them. The exphcit table of uncertainties (given as fractional errors) made the SKM work 
very useful, as these that could be readily input into BBN codes. 

An important recent study takes a very different approach to the BBN error budget. Buries, 

Nollctt, Truran, & Turner (1999) and Nollett & Buries (2000) made an extensive Monte Carlo 
calculation of the light clement abundances and their uncertainties. In this work, the R- and S- 
factors for the 11 key strong reactions are generated from fits to simulated data sets. These mock 
data are drawn randomly, on the basis of the actual data. For each realization of the nuclear data, 
the thermonuclear rates are computed via eq. (1), and run in the BBN code. Care is taken to 
allow for correlations due to normalization errors. We note however, that there is a danger in this 
procedure, as the resulting adopted S and R- factors may not have a dependence as a function of 
energy as expected from theory. Indeed, some cases, such as ^Li(p,a)*^IIe, show a behavior which 
is probably not physical and is simply a result of the available data at particular energies^ 

Thus, Nollett & Buries do not directly report a mean and error for the thermonuclear rates 
(though they do produce mean values and errors for the R- and S'-factors). Instead, they report 
the final results for the abundance mean values and errors, as a function of f]. Finally, we note that 
Nollett & Buries adopt a neutron lifetime = 885.4 it 2 s, which is slightly lower than the current 
world average of the Particle Data Group (Groom et al. 2000), 885. 7 it 0.8 s. This difference affects 
mostly ^He, with a small drop AYp c± —0.0001 for the Buries Sz Nollett value, although the error 
used here is considerably smaller. 

To compare our nuclear error budgets with those of SKM and Nollett & Buries (2000), we have 



^This was pointed out by A. Coc in his talk at Cosmic Evolution, Paris, Nov. 2000 (Coc, Vangioni-Flam & Casse 
2001). 
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plotted fractional errors in the R- and S'-factors for 10 of the 11 strong reactions.^ These appear 
in Figure 3. We see that our sample variance uncertainties are always tighter than those of SKM. 
On the other hand, our results are almost always larger than or about equal to the errors found by 
Nollett &: Buries.^ It is interesting that, with the exception of the two (0,7) reactions, the Nollett 
&; Buries fractional errors are roughly constant with energy. In four of these cases, our values 
agree well with the average of the range spanned by Nollett &; Buries. For 7)^He, our value 
is slightly higher, due to method's slightly more conservative treatment of the discrepant data for 
this reaction as is the case for '^Li{p, Q!)'^Hc. In the mirror reactions of t{d, n)^He and '^He((i, p)'^He, 
the reactions proceed through a resonance, and thus the errors are very sensitive to the goodness 
of fit for the NACRE and SKM curves, which is not an issue for the empirical approach of Nollett 
&; Buries. 



4. The Impact on BBN 

We have now adopted a set of thermonuclear rates using the NACRE compilation (§2), and 
we have presented several estimates of the uncertainties in theses rates (§3). In this section, we will 
put these rates into the BBN code, and compiitc the light element abundances and their errors as 
a function of r]. We can then quantify the impact of the NACRE compilation on the mean values 
of the abundances, and on their uncertainties. 

The implementation of the Monte Carlo is standard and well-described in, e.g., Krauss & 
Romanelli (1990) and SKM. Briefly, a grid of t] is chosen, with values in the (1 — 10) x 10^^*^ range. 
At each ij, the code is run 1000 times. For each run, we generate 12 Gaussian random numbers Zi 
with zero mean and unit variance; these are used to choose the thermonuclear reaction rates Aj in 
terms of their mean values Aj^o a^nd fractional errors fi = SXi/Xifi-. 

\i{T) = Xi,o{T) [1 + ZiMT)] (15) 

Where fi{T) is constant for the minimal uncertainty and sample variance cases, and T-dependent 
for the high/low case. 

At each rj, the central value of the light element abundance yi is taken to be the mean of the 
1000 values from the runs, and the standard deviation a{yi) is taken from the sample variance of the 
1000 runs. A technical note: the sample variance itself has an expected deviation of Sa"^ ~ a'^/N 
from the true error, so that our uncertainty estimates have a statistical error of l/\/iV ~ 3%. To 
allow for this, we increase the computed sample variance by 3%; in every case the effect is small. 



®For legibility, we have omitted p{n, ■y)d, for which Nollett & Buries adopt a constant fractional error of 10% at 
the 95% CL. The corresponding SKM uncertainty is 14% and our sample variance error is ~ 9% (see Table 2). 

'^Note that the plotted Nollett & Buries errors are 1/2 of their 95% CL uncertainties. 
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4.1. NACRE Central Values 

It is important to first compare the direct impact of the NACRE rate compilation on BBN 
predictions. To do this we will use the expanded BBN code of Thomas et al. (1993), which was 
based on the Kernan (1993) code, with the Monte Carlo implemented by Hata et al. (1996). The 
uncertainties in the 1996 Monte Carlo were derived from SKM. The central values of the light 
element abundances appear in Figure 4a, which shows results for the NACRE rates along with 
the Hata et al. (1996) results, (the solid and dashed curves). We see that the two predictions arc 
almost identical to each other. The differences between the two cases are plotted in Figure 4b as a 
percentage of the 1996 predictions. We see that the changes are small: the largest shift is a 10% 
rise in D/H at 77 = 10~^, and that below rj = 3x 10""^'^, all shifts are smaller than 5% in magnitude. 

These results confirm the findings of Vangioni-Flam, Coc, & Casse (2000): the central values 
of the BBN predictions are insensitive to the reaction network choice. It is both non-trivial and 
comforting that different rate tabulations, with different functional forms, yield the same central 
values. This agreement reaffirms that basic predictions of BBN are robust. 

Although the central light element predictions of the two rate compilations do not differ sig- 
nificantly, a larger effect is found by renormalizing the rates to reflect the data in the BBN energy 
range. Figure 5a shows the predictions when the rates are renormalized as described in the previous 
section. We see that central values do change, at levels larger than those discussed above. Figure 
5b shows the percent difference of these predictions with those of Hata et al. (1996). We see that 
the shifts are ^ 1% for Yp, at a level of < 15% for D and ^He, and vary the most for ^Li, which 
runs between ±20%. The shifts trace back to the rates which have the largest renormalizations. 
In particular, the d{p, 7)'^He rate drops by 7.6% with respect to NACRE due to the preferential 
weighting of the Schmid ct al. (1995) data with their smaller errors (Figure Id). This change alone 
is predominantly responsible for the shifts in D, ^He, and ''Li at 77 > 3 x lO"^'^; the ^He((i,p)^He 
-|-6.9% renormalization is responsible for most of the rest of the shift in this regime. At smaller rj, 
there is a rise in ^Li; this is predominantly due to the combination of the increase in ^Li production 
via a 5.8% rise in the t{a, ^Yhi rate, and a decrease in '''Li destruction due to a 4.5% lowering of 
the ''Li(p, a)^He rate. We note that these shifts in the central values of the predicted abundances 
are well within the calculated uncertainties discussed below. 



4.2. NACRE High/Low Errors 

The impact of NACRE's high/low errors (§3.2) on the light element abundances is seen in 
Figure 6a. The logarithmic nature of the standard abundance plot, and the rapid variation of the 
trace elements, makes it difficult to assess and compare the errors. Thus we have also plotted the 
fractional errors ai/ni in Figure 6b. We see that overall, the NACRE high/low errors are quite 
comparable to those found using the SKM estimates. The most significant effect is the reduction in 
the ^Li errors, which is due to the lower uncertainties in the mass-7 production reactions i(a,7)^Li 
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and ^He(a,7)^Be. Here, the errors are a strong function of 77 owing to the ^Li and ^Be production 
channels, and the errors are comparable to those of SKM at 7710 = lO^^rj near 3. However, the error 
on Li is reduced by as much as ~ 25% for rjio near 1 or 10, 

We confirm here the key result that the theoretical errors in BBN are still significant. The 
errors in D and ^Li in particular are only marginally smaller than the accuracy of current or near- 
future observations. Thus, improvements in the nuclear inputs to BBN can still have an important 
effect on cosmology, as we will discuss below (§5.3). 

It is important to note that our results differ somewhat from those of Vangioni-Flam, Coc, 
& Casse (2000), due to differences in approach. Vangioni-Flam et al. used the high/low errors to 
examine the effect on BBN of each reaction individually. They then found the abundances when 
all the reactions are set to their high limits, then when all reactions are set to their low limits. This 
procedure is a simple and rapid means of gauging the net effect of the uncertainties. However, as 
Vangioni-Flam ct al. note, this can lead to compensating effects (e.g., since both the ^Li destruction 
and production rates are raised, there is a degree of cancellation in the effect of the errors on the 
final abundances). Thus, the Vangioni-Flam et al. net uncertainties are lower limits to the true 
uncertainties; this fact is borne out in Figures 6a and 6b, where the errors are equal to or larger 
than those of Vangioni-Flam et al. 

4.3. Minimal Uncertainties: Ax^ = 1 

We now turn to the case of the uncertainties given in the Ax^ = 1 analysis of §3.3. Results 
appear in Figure 7, where we see that the errors are now very small. This is not a surprise, given 
that these errors represent a lower limit to the true uncertainties. In this limit, the nuclear input 
errors are small enough to be negligible in comparison to the uncertainties one can expect in the 
observed abundances. 

The errors presented here are derived by assuming that all of the cross section data points 
may be combined without concern for systematic errors between the experiments. As this is not 
likely to be the case, the true BBN theoretical error budget is considerably larger, however. We 
now turn to the "sample variance" error estimator, which allows for these effects. 

4.4. Sample Variance 

The BBN predictions based on the sample variance errors appear in Figure 8. The smallest 
errors are of course for ^He, where the uncertainties are of order ~ 2%, or Ay = 0.0005. We 
note that the neutron lifetime is now sufficiently well-known that the prediction for '^He are now 
dominated instead by the nuclear errors. The D/H errors are also of the same order as the statistical 
errors in the (high-redshift) data, though systematics are likely to dominate the observational 
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uncertainties, as discussed in the next section. 

The largest errors are for ^Li, where the ~ 10 — 18% uncertainty is of the same order of 
the statistical error in the observations. Again, systematic errors are more important for the 
observations. Even so, it is particularly important to reduce the error in Li regardless of the 
uncertainty in the observations. As noted by, e.g., Ryan et al. (2000), since the interesting range 
in r) places the Li(?7) curve near its minimum, it is therefore slowly varying. Even if there are zero 
observational errors, the shallowness of the ^Li — r) curve means that the theoretical uncertainty in 
the Li prediction leads to a large error in the allowed range for 77 (below, §5). 

We now compare the errors derived from the sample variance with the other estimators we have 
used. Referring to the panels (b) of Figures 8, 6, and 7, we see that the sample variance errors are 
smaller than those using SKM, or NACRE's high/low ranges, but larger than the minimal errors, 

as expected. Relative to SKM, the Li errors are reduced on the low-r/ side (i.e., the ^Li regime) by 
the inclusion of the precision t(Q;, 7)''Li, data of Brunc, Kavanagh, & Rolfs (1994). On the high-r/ 
side, the reduced Li errors follow from our smaller errors in '^He(a, 7)''Be. The '^Hc and D errors 
are considerably lower (by roughly a factor ~ 2), due to the smaller errors in p{n,'y)d, d{d,n)^ile, 
and d{d,p)t. 

Of the error estimation methods we have used, we feel the sample variance method provides 
a good and fair estimate of BEN uncertainties. This procedure has the virtue of allowing for 
systematic uncertainties due to reaction normalization errors. It also has the virtue of simplicity 
of definition and implementation, and one can judge by eye the goodness of the fit. 

The sample variance errors are compared with those of Nollett Sz Buries (2000) and SKM in 
Figure 9. The differences among the fractional errors plotted stem mostly from differences in the 
errors ai themselves, except for the case of ''Li, where the differences in central values (ours are 
lower than those of SKM and NB) also contribute to changes in a/fi. We sec in Figure 9 that 
our relative errors generally lie between those of SKM and Nollett & Buries. This is as expected, 
since the NACRE thermonuclear rates, normalized or not, follow similar trends to those rates used 
by SKM and Nollett & Buries. Also, our sample variance usually lies between those of SKM and 
Nollett &; Buries' errors. For each element, different reactions dominate the errors at different 
regimes in 77, and thus the variations of the shapes of the curves in each panel of Figure 9 reflect 
the differences in the underlying reaction errors. As seen in Figure 3 and Table 3, both the SKM 
and Nollett & Buries errors are generally of comparable magnitudes across most reactions; ours 
have somewhat larger variations in magnitudes, and thus our curves have a somewhat different rj 
dependence than those of these compilations. 

For example, in comparison with Nollett &; Buries, we find that the D/H fractional errors are 
in very good agreement for rjio < 3. At higher rj, our errors then grow larger, about 30% larger 

at ??io = 5, to 80% larger at ryio = 10; this difference is mostly due to the increasing sensitivity of 
deuterium to the d{p, 7)^Hc rate, for which our adopted error is larger. The trend for ^He shows 
a more pronounced rise at higher r], due to the increasing importance of both the d{p, 7)^He and 
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^He(c?,p)^He rates, where again we have used a larger uncertainty. 

For ^Li, our errors are consistently higher, by a factor that varies from 70% to 100%; part of 
this difference arises from the somewhat higher central values Nollett & Buries have at higher rj 
(and thus a smaller cr/ii). The rest of the differences trace back to the systematically tighter Nollett 
& Buries error bands, as discussed above. The difference docs not stem from a single reaction, but 
rather the cumulative effect of the lower Nollett Sz Buries errors for several reactions: d{p, 7)'^He, 
as well as the mass-7 producing and destroying reactions. In the case of Y, our errors are smaller 
than those of the other groups; this arises because we have adopted the new errors in the neutron 
mean life which are smaller than those used by Nollett & Buries and others in the past. 

5. Comparison with Observation 

The procedure for comparing the light element abundances as predicted by BBN theory with 
the light elements as observed astrophysically is well-described elsewhere, e.g., (Olive, Steigman, 
& Walker (2000)). Here, we simply summarize the needed observational inputs and statistical 
techniques. 

5.1. Observational Inputs 

Data for ^He is available from about 70 extragalactic H II regions (Pagel et al. 1992; Skillman 
& Kennicutt 1993; SkiUman et al. 1994; Izotov, Thuan, & Lipovetsky 1994,1997; Izotov & Thuan 
1998b) and has been recently compiled in Olive, Steigman, k, Skillman (1997) and Fields & Olive 
(1998). Unfortunately, the determination of the ^He abundance in these systems is not particularly 

straightforward. To convert the observed He emission line strengths into abundances requires 
knowledge of several physical parameters describing the H II system, such as the temperature, 
electron density, optical depth, and degree of underlying stellar absorption. 

The older data of Pagel et al. (1992), Skillman &: Kennicutt (1993), and Skillman et al. (1994), 
used S II data to fix the electron density from which was derived a relatively low ^He abundance. 
Higher ^He abundances were derived by Izotov, Thuan, k. Lipovetsky (1994,1997) and Izotov k. 
Thuan (1998b), where the electron density was determined self consistently from five distinct He 
emission lines. In addition, it was pointed out in Izotov k Thuan (1998a) that one of lowest 
metallicity regions, I Zw 18 NW, is plagued by underlying stellar absorption. The compilation 
by Izotov k Thuan (1998b) of 45 extragalactic H II regions yielded a primordial ^He abundance 
(by performing a regression on the data versus the 0/H abundance) oi Yp = 0.244 zt 0.002. As 
the compilation of Izotov k Thuan (1998b) also presented results based on S II densities, it was 
possible to combine all of the data in a systematic way. The result of Fields k Olive (1998) 
which included the data of Pagel et al. (1992), Skillman k Kennicutt (1993). and Skillman et 
al. (1994) and Izotov k Thuan (1998b) (but without the possibly erroneous I Zw 18 NW), yielded 
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the somewhat lower result of Yp = 0.238 ± 0.002. Note that the data of Izotov & Thuan (1998b) 

alone based on S II densities give 0.239ib 0.002. One should also note that a recent determination by 
Peimbert, Peimbert, &: Ruiz (2000) of the ^He abundance in a single object (the SMC) also using 
a self consistent method gives a primordial abundance of 0.234 it 0.003 (actually, they observe 
Y = 0.240 ± 0.002 at the relatively high value of [0/H] = -0.8, where [0/H] refers to the log of the 
Oxygen abundance relative to the solar value). 

Recently, a detailed examination of the systematic uncertainties in the ^He abundance deter- 
mination was made by Olive &; Skillman (2000) . There a Monte Carlo simulation of synthetic data 
showed that the He abundance determinations using a minimization in the self-consistent method 
typically under-estimated the true errors by about a factor of 2. The reason for the enhanced error 
determinations is a degeneracy among the physical parameters (electron density, optical depth, and 
underlying stellar absorption) which yield equivalent results. Indeed, it was found that there 
are competing biases in the data. On the one hand, the presence of underlying stellar absorption 
(expected to be present to some level in most systems) leads to an underestimate of the He abun- 
dance. While on the other hand, the minimization solutions tend to find solutions with erroneously 
low densities thus minimizing the collisional corrections and over-estimating the He abundance. 

With these cautions in mind, we will never-the-less adopt the Fields & Olive (1998) value of 

Yp = 0.238 ± 0.002 ± 0.005 (16) 

where the first error is statistical and the second systematic. Clearly, the dominant error source is 
systematic, and while we have given our best estimate for its value, this number may yet turn out 
to be too small until a complete reanalysis of the existing or new data is possible. For the full error 
in Yp, we will use the quadrature sum of the two components. 

As in the case of ^He, there is a considerable body of data on ''Li. There are well over 100 hot 
halo dwarf stars with ^Li observations. The discovery of a plateau in the Li abundance for hot and 
metal poor dwarfs by Spite Sz Spite (1982), is generally recognized as the primordial value. Recent 
high precision studies of Li abundances in halo stars have been made by Bonifacio & Molaro (1997) 
and Ryan, Norris, & Beers (1999) which confirm the plateau. The latter in fact identified a small 
but significant Li-Fe trend in the low metallicity regime. Ryan et al. (2000) showed that this trend 
is indeed expected given the presence of cosmic ray interactions in the early Galaxy. These provide 
a source of Galactic Li production on top of a cosmological value (very analogous to the derivation 
of Yp from a regression ofY — Z data and the dY/dZ slope). Ryan et al. inferred a primordial Li 
abundance of 

^Li/H = (1.23l°:|) X 10"^° (17) 

Note that correcting for Galactic production lowers ^Li/H compared to taking the mean value over 

a range of metallicity. For our adopted primordial ''Li, we assume a gaussian distribution with 
central value equal to that of eq. (17), and a standard deviation which we conservatively take to 
be the larger of the asymmetric errors. 
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We note that in contrast to the downward correction due to post big bang production of Li, 
there is a potential for an upward correction due to depletion. While most studies of possible 
depletion factors are small, they may not be negligible. Both Vauclair & Charbonncl (1998) and 
Pinsonneault et al. (1998) argue for an upward correction of 0.2 dex. We note however, that 
the data do not show any dispersion (beyond that expected by observational uncertainty) and 
both Bonifacio & Molaro (1997) and Ryan, Norris, & Beers (1999) have argued strongly against 
depletion. We further point out that the observation of the fragile isotope, ^Li, (Smith, Lambert 
& Nissen 1992, 1998; Hobbs and Thorburn 1994, 1997, Cayrel et al. 1999; Nissen et al. 2000) also 
puts strong constraints on the degree of depletion (Steigman et al. 1993; Lemoine et al. 1997; Fields 
& Olive 1999; Vangioni-Flam et al. 1999). 

The observational status of primordial D is complicated but also promising. Deuterium has 
been detected in several high-redshift quasar absorption line systems. It is expected that these 
systems still retain their original, primordial deuterium, unaffected by any significant stellar nucle- 
osynthesis. At present, however, there are four good determinations of D/H in absorption systems 
with considerable scatter. Two good measurements (Buries & Tytlcr (1998a; 1998b)) with D/H 
determinations of D/H = (3.3 ± 0.3) x 10"^ and (4.0 ± 0.7) x 10"^ have recently been supple- 
mented by the new results of O'Mcara et al. (2000) with D/H = (2.5 ±0.2) x 10"^ These authors 
measured D/H in a Lyman limit system which has a neutral hydrogen column density a factor of 
30-100 higher than the previous D/H systems, and in which most of the H is neutral rather than 
ionized. As a result, O'Meara et al. (2000) were able to detect D in five transitions. The inferred 
D/H value is lower than the previous absorption systems, and the three systems combined give 

D/H = (3.0 ± 0.4) x 10-^ (18) 

O'Meara et al. (2000) note, however, that = 7.1 for the combined 3 D/H measurements (i.e., 
v = 2), and interpret this as a likely indication that the errors have been underestimated. 

Significantly more discrepant is the higher determination of D/H, (2.0 it 0.5) x 10~^, as seen 
in one low-redshift system (Webb et al. (1997); Tytler et al. (1999)). Clearly there are systematic 
uncertainties which are yet to be sorted out. Indeed, Levshakov, Kegel &: Takahara (1998) have used 
the data in the high redshift system of Buries &; Tytler (1998a) but with a different model for the 
velocity distribution of the absorbing gas, and derived a 95% confidence range 3.5 xlO~^ < D/H < 
5.2 X 10~^. Levshakov, Tytler & Buries (1998) have also applied this model to a reanalysis of the 
system of Buries & Tytler (1998b), finding a a 68% confidence range of D/H ~ (3.5 — 5.0) x 10~^. 
Because of these complications, we will examine the impact of low and high D on BBN. 

5.2. Likelihood Analysis 

To quantitatively compare the observed abundance data with theory, we use a likelihood anal- 
ysis. The formalism is described elsewhere (Fields &; Olive (1996); Fields, Olive, Kainulainen, k, 
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Thomas (1996)), but the basic idea is that the Monte Carlo results (means, variances, correla- 
tion coefficients) allow one to compute the theory likelihoods for the elements as a function of rj. 
The convolution of these with the observed abundance distributions gives the overall likelihood 
as function of r/; these can determine the goodness of fit of the model and allow for parameter 
estimation. 

Figure 10 shows the likelihoods derived using the sample variance errors. In panel (a), we see 
the theory-observation combined distributions for individual elements. The peaked nature of the 
^He and D curves stems from the monotonic nature of the 77 dependence of these abundances. The 
flatness of the ^Li likelihood arises because the observed ^Li abundance we have adopted lies at the 
minimum of the lithium trough, so that the usual two-peaked structure has merged into a single 
broad feature. As has noted in Fields et al. (1996) and elsewhere, the ^He and ^Li likelihoods 
are in excellent agreement with each other and with the high D distribution. This is reflected 
in the combined likelihood distribution (panel b). Here, L47 represents the combined likelihood 
distribution using only ^He and ^Li, while L247 is the combined likelihood distribution using D, 
^He and ^Li. A value of order unity is an indication that the agreement is good. 

On the other hand, the low D distribution agrees marginally at best with the ^He and ^Li data. 
This discrepancy has been widely noted, and if real, suggests that either one or more of the observed 
abundances is incorrect, that the systematic uncertainties have been grossly underestimated, or 
perhaps that new physics is at play. One should recall, though, that the D observations could 
as well suffer from systematic errors that are not reflected in numbers we have adopted. Indeed, 
systematic uncertainties are not included in the error budget of eq. (18). Since D/H is a strong 
function of rj, the D likelihood distribution is very sensitive to changes in the D data; for example, 
a modest shift upward in the D central value, or a larger D error, could lead to agreement with 
^He and ''Li. For example, even if D/H is only as high as ~ 5 x 10~^, as deemed possible in the 
mesoturbulent models of Levshakov, Kegel & Takahara (1998) and Levshakov, Tytler & Buries 
(1998), the peak of the likelihood function shifts down to 7710 ~ 4 and would indeed be in good 
agreement with the other light elements. As is also well known, the high D/H distribution is 
in excellent agreement as seen in Fig. 10. Clearly, a determination of the true, primordial D 
abundance-and an accurate assessment of its errors — is of paramount importance. 

Had we used the high/low errors of the NACRE rates we would have obtained very similar 
likelihood distributions and predictions for 77. Similarly, the use of the unrenormalized rates would 
not have made a big difference here cither. In the extreme case that wc use our minimal uncertainties 
based on Ax^ = 1, we would of course have tighter predictions as shown in Figure 11. It is 
interesting to note that even in this case, the likelihood distributions for "^He and ^Li are broad due 
to residual systematic uncertainties in the observations. Therefore, for any significant improvement 
in the BBN predictions, the systematic uncertainties will need to be resolved. 
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5.3. Constraints on Cosmology: The Baryon Density and Light Neutrinos 

Using the results of the hkehhood analyses in the previous section, we can derive limits on rj. 
We will focus on the sample variance case; as we have remarked, the likelihoods for the high/low 
case are only slightly different, and we have verified that they give very similar results for rj. 

For the sample variance error, we derive a most-likely value of ry, i), and a 95% CL range. For 
the ^He-^Li combined likelihood, we find (to two significant digits) 

r)io = 2.4 Ob/i^ = 0.0089 (19) 
1-7 < mo < 4.7 ^ 0.006 < Ob/i^ < 0.017 (95% CL) . (20) 

where rj is the value at maximum likelihood. We use the conversion r/io = 274i7B^^- For high D, 
the result is very similar, owing to the close overlap of the individual element likelihoods; we have 

= 1.9 ^ n^h^ = 0.007 (21) 
1-6 < 7710 < 3.3 0.006 < Ob/j^ < 0.012 (95% CL) . (22) 

For low D, the goodness of fit is worse, as we have noted, and we have 

fjio = 5.3 ClBh^ = 0.019 (23) 
4.7 < r?io < 6.3 0.017 < ^Bh^ < 0.023 (95% CL) . (24) 

We have of course used the standard model in the above determinations of rj. In particular, 
we have assumed N,y = 3. As was shown in previous analyses (Olive &: Thomas 1999, Lisi, Sarkar, 
& Villante 1999), the combination of ^He and ^Li (as well as high D/H) is quite compatible with 
the standard model value, and we do not expect that result to change here. Similarly, the result 
including low D/H is at best compatible with ATj, = 3 at the 2 a level. A complete two-dimensional 
likelihood analysis based on the NACRE compilation is beyond the scope of the present work. 



6. Constraints from Cosmology: The Microwave Background 

BEN has long provided the best estimate of 77, and thus the cosmic baryon content (i.e., 

Jl^/i^). However, measurements of the anisotropy spectrum of the CMB can constrain many key 
cosmological parameters, including 77 (see e.g.. White, Scott, & Silk (1994)). Thus, the agreement 
between the two estimates of cosmic baryons provides a fundamental test of cosmology (Schramm 
& Turner (1998)). 

Recent anisotropy data already allows for an initial investigation of the range in rj favored by 
CMB. The acoustic peaks in the CMB at small angular scales are sensitive to the baryon density 
as well as to other parameters, and the most reliable estimates of 77 come from fitting over a range 
of angular scales (including multiple peaks and valleys). The BOOMERanG-98 (de Bernardis et 
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al. (2000)) and MAXIMA-1 (Hanany et al. (2000); Balbi et al. (2000)) experiments provide the 
first data wcU-suitcd for tliis effort, as they span the entire first acoustic peak and the region where 
the second peak is expected. 

For these experiments, the inferred value of r/ is sensitive mostly to ratio of first to second 
peaks. The second peak appears weakly, if at all, in the BOOMERanG-98 data, and somewhat 
more strongly in the MAXIMA-1 data. The two groups combined their data (and those of COBE), 
and made estimates for a set of cosmological parameters (Jaffe et al. (2000)). They estimate 
^Bh? = 0.032^0^1 (68% CL), which corresponds to 

7?CMB = (8.8li:t) X 10-1° (68% CL) (25) 

Similar analyses by other authors give similar results^ (e-g-, Tegmark & Zaldarriaga (2000); Tegmark, 
Zaldarriaga, & Hamilton (2000)). 

The CMB preferred range in rj (eq. 25) is about 3 a away from the BBN results derived from ^He 
and ^Li (Figure 10b and eq. 19). Including low D/H only slightly reduces the difference, to two a. 
Indeed, the CMB favored range based on the BOOMERanG-98 and MAXIMA-1 data does not agree 
with any of the three light nuclides. For example, the CMB data demand that D/H = 1.6 x 10"^, 
which is significantly below all high-redshift determination of deuterium. Moreover, the CMB- 
preferred D/H lies just at the present-day local interstellar value, D/H = (1.5 it 0.1) x 10"^ (e.g., 
Linsky et al. (1998); Sahu et al. (1999)). This would allow for almost no processing of D/H over 
the history of the Galaxy. This is very unlikely even if primordial infall played a big role in the 
chemical evolution of our Galaxy. The CMB data also require Yp = 0.252 and ''Li/H = 7.4 x 10^''. 
For ^He, this is somewhat high, but given the magnitude of the systematic uncertainties, it would 
be difficult to exclude this value outright. On the other hand, the ^Li value would require a factor 
of ~ 6 in depletion. This is probably a factor of ~ 3 too high relative to the most optimistic models 
of stellar depletion. 

If the discrepancy is real, it may point to new physics occurring sometime prior to recombi- 
nation. Kaplinghat & Turner (2001) pointed out that entropy production after BBN but before 
recombination could account for the difference without requiring new physics at the BBN epoch. 
On the other hand, new BBN physics could also lead to a higher baryon density for the same 
abundance constraints. In particular, there has been considerable attention recently to the case of 
BBN with a large neutrino chemical potential (Esposito, Mangano, Melchiorri, Miele, & Pisanti 
(2001); Lesgourgues & Peloso (2000); Orito, Kajino, Mathews, & Boyd (2000); Esposito, Mangano, 
Miele, & Pisanti (2000b)). Finally, Kneller, Scherrer, Steigman, k Walker (2001), show that while 
the BBN-CMB discrepancy is robust given the BOOMERanG-98 and MAXIMA-1 data, the non- 
standard parameters depend sensitively on CMB priors. 

However, before we abandon the standard model, we note also CMB results of the two exper- 



In the cases where nucleosynthesis information was not assumed as a prior. 
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iments independently have very different second peak points, leading to different (but consistent) 
baryon contents. Furtliermore, tfie recent data from the Cosmic Background Imager (CBI; Padin 
et al. (2000)) imphes a very different result. CBI is sensitive to the higher order acoustic peaks 
{i = 400 — 1500), and their first results favor O^/i^ = 0.009, or r/io = 2.46. While this is in sharp 
contrast with the results of BOOMERanG and MAXIMA, this value is lower than the ij range 
suggested by the low D/H data. In fact, the CBI result is in excellent agreement with the rj range 
favored by ^He and ^Li (Figs 10). Clearly, more and better CMB data are needed to clarify this 
situation. 

Fortunately, high-quality CMB data will soon be at hand. The MAP satellite (Wright (1999)) is 
scheduled to launch in 2001. and will return data about a year later. These precision measurements 
should give to better than 10%. The PLANCK Explorer, scheduled for launch in 2007, will 
reduce these errors to better than 3%. As these data become available, the comparison of the rj 
predictions with those of BBN will become an acute test of cosmology. Clearly, a strong discrepancy 
will be a surprise and will demand an explanation. 

On the other hand, it is possible that, with precision results available for both BBN and the 
CMB, the baryonic predictions will be in agreement. Then, as the CMB results improve, it will 
become useful to use the CMB range for as an input to the BBN analysis, and to use this, e.g., 
to give precision predictions for the light element abundances. An illustration of such a prediction 
is shown in Figure 12, which shows the light element abundance predictions assuming a constant 
central value of rjio = fjio (as in eq. 19), but with an error which goes from that of our current 
limits down to 10% (MAP) and then 3% (PLANCK). We see that these experiments will be able 
to determine the primordial abundances quite accurately and thus will impact studies of galactic, 
stellar, and cosmic-ray evolution. 

A noteworthy feature of Figure 12 is the slight shift in the peak values of the likelihoods as the 
limits on rj become tighter. Thus, the maximum likelihood estimator for the abundance y is slightly 
offset from yBBN('?) but converges to this value as the CMB uncertainty in rj becomes small. This 
comes about because the plotted likelihood is given by the convolution 
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The BBN theory hkelihood 



Lbbn oc cJbbn exp[-(i/ - yBBN)^/2cr|BN] 



(27) 



depends on r/ in a complicated way through the abundance and error curves ?/bbn(??) and o"bbn('?)- 
One can show that the peak of L{y) is shifted due to the variation of the mean values and errors 
with T]. One can also show that the shift goes to zero as the CMB error becomes small. 



-26- 



7. Conclusions 

We have examined the effect of the NACRE thermonuclear reaction rates on BBN. We verify 
that the central values of the abundance predictions are very similar to those used in previous 
studies (cf. Hata et al. 1996). We note that the NACRE fits are in general good representations 
of the shape of the data, but the normalizations do not precisely minimize x^- We have thus 
computed the renormalizations needed to minimize x^- Although small, the renormalizations do 
lead to changes with the previous predictions, particularly in that the Li prediction at high ij is 
reduced, and thus less Li depletion is required in this regime. 

We have used the NACRE rates and database to determine the theoretical uncertainty in 
BBN via the standard Monte Carlo procedure. To do this we have examined the effect of different 
estimates of the nuclear reaction uncertainties. NACRE's high/low limits give BBN uncertainties 
similar to those of SKM. Using a simple sample variance procedure, we arrive at errors which are 
smaller than these, while being slightly larger than those of Nollett and Buries. We have also 
computed lower limits to the errors on the reaction rates, which lead to the minimal theoretical 
errors on the light element abundances. 

Error studies such as ours are increasingly important as cosmology enters the precision era. In 
particular, an independent prediction of rj is emerging from CMB anisotropy data. At present, the 
CMB results are tantalizing, and possibly discrepant with the BBN predictions, though the CMB 
predictions for r] seem to be uncertain at the moment. At any rate, the CMB data will improve 
dramatically with the upcoming launch of MAP; this experiment, along with PLANCK, will make 
a strong test of BBN and of cosmology. 

In anticipation of this test, further work is needed to find a realistic error budget for both 
the theory and observation of light element abundances. A key input on the theory side can 
come from accurate nuclear experimental data and carcfiil analyses of the error budget in key 
reactions. Specifically, we urge more data be taken for p(n, 7)^, in order to confirm theory, and put 
this reaction on a solid empirical basis. Other key reactions are those with internally discrepant 
reactions, d{p,^)^Ile and ^He((i,p)''He, and the lithium sources and sinks, i(a,7)^Li, ^He(a,7)^Be, 
and ^Li(p, Q;)^He. 

We are grateful to the NACRE collaboration, and particularly Carmen Angulo and Pierre 
Descouvrement, for their helpful responses to our questions about the data and rate compilations. 

We are indebted as well to A-Iichacl Smith for instructive discussions and sharing reaction fitting 
functions. We thank Sam Austin, Carl Brune, Michel Casse, Alain Coc, Vijay Pandaripande, and 
Elisabeth Vangioni-Flam for instructive conversations. 
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Fig. 1. — Experimental data and R- and S'-factor fits for four key reactions. The solid curves are 
our renormalized fits, with our sample variance error bands given by the short dashed curves. For 
the S-factor curves, data are given in NACRE, whose raw fits are given by the long dashed curves. 
a) The i2-factor for ^He(n,|?)t. The data are those of SKM with the addition of Brune et al. (1999). 
The fit is ours, as described in §2.2. b) The S'-factor data for 7)^He. c) The S'-factor data for 
^He(Q,7)^Be. d) The 5-factor data for t(a,^yLi. 
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Fig. 2. — NACRE "high/low" fractional errors for the seven NACRE reactions important for BBN. 
Dotted curves: high/low errors; solid curves: sample variance errors; dashed curves: SKM errors 
which are shown for comparison. 
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Fig. 3. — Fractional Icr errors in the R- and S-factors for sample variance analysis, as well as 
Nollett k. Buries (2000). Solid curves: sample variance limits; dashed curves: Nollett & Buries 
upper limits; dotted curves: Nollett & Buries lower limits. The p(n, 7)^ reaction is not shown, as 
Nollett & Buries adopt a constant la fractional error of 5% in this case. 
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Fig. 4. — Light element abundances as a function of rj, using NACRE rates without renormalization. 
a) The sohd curves are the central values of the NACRE predictions, dashed curves are those of 
SKM. The two are extremely close, so that they often appear to overlap, b) Percent difference 
100(yf^^c^^/y|996 _ 1) between the NACRE and Hata et al. (1996) central values. 
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Fig. 5. — As in Figure 5, for light element predictions with renormalization. 
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Fig. 6. — Light element abundance predictions, and their uncertainties, as a function of rj. a) The 
soHd curves are the (renormahzed) NACRE abundance predictions, and the broken curves are the 
la Monte Carlo errors for the NACRE high/low error estimates, b) The fractional errors in the 
light element predictions, for the NACRE high/low estimates as in (a) (solid curves). The fractional 
errors plotted are crj/^Uj, where fii is the mean value of abundance i, and ai is its error. 




Fig. 7. — As in Figure 6, with errors for the Ax^ = 1, "minimal uncertainty" estimates. Note the 
change in vertical scale between panel (b) and Figure 6b. 




Fig. 8. — As in Figure 6, with sample variance error estimates. 



- 39 - 




Fig. 9. — A comparison of BBN fractional errors by different groups. Solid curves show our sample 
variance errors for the four light elements, short dashed curves are the errors from Nollett & Buries 
(2000), and long dashed curves are the errors from SKM. We see that our errors generally fall 
between these two cases. 
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Fig. 10. — Likelihoods for the sample variance errors, a) Curves for ^He and ^Li, as well as high and 
low D individually. Plotted are the likelihoods for combined theory and observation distribution. 
The small amplitude for the low D case indicates a poorer fit. b) Total 2- and 3-element combined 
likelihoods. 
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11. — Likelihoods as in Figure 10 for the minimal errors. 



- 42 - 





'.225 0.24 0.255 
Y 



0.5 1 1.5 2 
lO^xD/H 



0.5 p 




0.3 

o 
o 

6 0.2 
. I — I 



o 







1 M M 1 


1 1 M M 1 M M M M 1 M M M 1 

7] error - 




current _ 




- - MAP J 


'/ 
I 


; Planck 


: i 

- * 

1 1 1 Ut"! 1 1 


ill i~TTTrrn-U-H-i-i-^iX 



0.8 1.6 2.4 3.2 
lO^xSHe/H 



12 3 4 
lOiOx^Li/H 



Fig. 12. — Light element abundances as might be predicted by r] inputs from future CMB experi- 
ments. We have assumed a central value of r/io = ??io (eq. 19). Solid curve: 10% errors in r/ (MAP); 
Broken curve: 3% errors in rj (PLANCK). We see that predictive power of the CMB limits on rj 
will become quite strong. 



